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SECTION  1 
INTRODUCTION 


Potential  engine  damage  caused  by  the  ingestion  of  dust-laden  air  is  a 
serious  consideration  for  the  operation  of  gas  turbine  engines  on  unimproved  runways  or 
in  dusty  environments.  Gas  turbine  engines  are  routinely  tested  for  the  effect  of 
ingestion  of  solid  particles  according  to  the  procedures  of  Military  Specification 
MIL-E-5007D.  However,  the  material  used  in  MIL-E-5007D  is  basically  crushed  quartz 
at  a  very  moderate  concentration.  This  dusty  environment  is  not  appropriate  to  test 
for  all  realistic  situations  of  potential  engine  damage.  The  reality  of  this  issue  has 
been  dramatically  emphasized  by  recent  (1,2)  incidents  in  which  gas-turbine-powered 
air  transports  have  attempted  to  traverse  volcanic  clouds.  The  material  associated 
with  the  volcanic  clouds  and  other  dusty  environments  is  more  of  an  earth-like  material, 
is  often  encountered  at  concentrations  well  in  excess  of  those  used  in  MIL-E-5007D 
tests,  and  has  a  particle  size  distribution  different  from  MIL-E-5007D. 

Several  different  mechanisms  can  be  active  in  altering  engine  performance 
during  or  after  exposure  to  a  dust-laden  environment.  For  most  modern  engines,  the 
operating  line  is  trimmed  as  close  to  the  surge  line  as  is  possible.  This  operating  line 
has  a  built-in  margin  for  known  degradation  terms  such  as  inlet  distortion,  transient 
gusts,  and  normal  component  wear.  Encountering  a  particle  laden  environment  of  the 
type  of  interest  here  falls  into  the  category  of  an  unanticipated  degradation  and  may 
be  manifested  in  one  or  more  of  the  following  ways:  (a)  glassification  of  hot-section 
components,  (b)  erosion  in  compressor  blading  and  rotor  paths,  (c)  blockage  of  cooling 
paths,  and  (d)  oil  system  or  bleed  air  supply  contamination. 

If  glassification  occurs,  it  will  likely  take  place  on  hot-section  components 
of  the  machine.  Whether  or  not  it  will  happen  depends  upon  many  parameters,  including 
engine  operating  point,  the  particular  engine  construction,  and  the  physical  and  chemical 
properties  of  the  dust  cloud  material.  For  the  particular  engines  discussed  in  this 
report,  glassification  did  not  occur  because  the  turbine  inlet  temperature  was  below 
that  required  for  equilibrium  melting  of  the  material  utilized. 

Contamination  of  the  bleed  air  systems  was  experienced  during  the 
measurements,  but  these  results  will  not  be  discussed  herein.  The  interested  reader  is 
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referred  to  other  project  reports  (10-12).  Severe  compressor  erosion  was  also  experienced 
by  all  of  the  engines  used  in  this  measurement  program  and  this  particular  aspect  of 
the  result  will  be  described.  The  problem  of  component  erosion  by  particulate-laden 
air  flows  has  been  the  subject  of  study  for  many  years  (3-9)  and  is  known  to  result  in 
the  loss  of  surge  margin. 

The  discussion  in  this  report  will  be  confined  to  the  effort  expended  in 
predicting  the  influence  of  erosion  on  the  engine  performance  and  in  reconciling  these 
predictions  with  the  experimental  data  obtained  from  two  TF33  turbofans  and  one  357 
turbojet  exposed  to  the  same  dust  environment.  One  of  the  TF33  engines  and  the  357 
engine  that  were  exposed  to  the  same  dust  mixture  were  completely  disassembled  so 
that  detailed  post-test  measurements  of  the  components  could  be  obtained.  As  a  result 
of  the  detailed  measurements  performed  on  the  deteriorated  engines  and  the  operational 
experience  obtained  with  them,  it  is  possible  to  evaluate  the  predictive  ability  of  a 
model  developed  in  the  following  sections  of  this  report. 
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SECTION  2 

BRIEF  REVIEW  OF  EXPERIMENTAL  APPARATUS  AND  OPERATION 


The  TF33  turbofan  and  the  357  turbojet  discussed  in  this  report  were 
operated  in  DNA's  Large  Engine  Research  Cell  (LERC).  The  LERC  is  a  test  cell  sized 
for  high  thrust  large  weight  flow  gas  turbine  engines  and  is  described  in  detail  in  (10). 

2.1  TEST  CELL  AND  INSTRUMENTATION. 

The  engines  were  mounted  in  the  LERC  on  an  A/M-37T-6C  test  stand. 
The  cell  is  an  air-cooled  noise  suppression  structure,  composed  of  primary  and  secondary 
inlets,  an  augmentor  tube,  a  perforated  exhaust  basket  and  an  exhaust  stack.  The 
secondary  inlet  carries  air  entrained  by  the  hot  exhaust  of  the  engine  to  the  exhaust 
basket  where  it  passes  through  perforations  which  break  up  the  audible  energy  of  the 
flow.  The  secondary  inlet  also  serves  to  both  cool  the  exhaust  and  to  relieve  the 
pumping  load  in  the  test  chamber.  It  is  relatively  easy  to  create  significant  cell 
depression  and  ram  effect  biasing  of  engine  thrust  in  a  facility  with  only  a  primary 
inlet.  However,  because  LERC  has  both  primary  and  secondary  air  intakes  it  sustained 
a  small  cell  depression  of  -.03  psig  during  operation  of  a  TF33  at  takeoff  rated  thrust. 

The  Control  Cab  is  part  of  the  A/M-37T-6C  test  stand  which  produced  a 
sound  proof  environment  for  the  engine  operations.  The  Cab  is  equipped  to  provide  all 
of  the  engine  related  readings  that  the  pilot  would  see  in  the  cockpit.  The  experiment 
was  intentionally  configured  in  this  way  in  order  to  make  the  results  as  directly  relevant 
to  a  pilot's  viewpoint  as  possible  (10). 

The  nearby  control  room  contains  readouts  of  cell  and  some  engine 
parameters.  The  parameters  which  most  directly  indicate  the  engine's  condition:  EGT, 
N2,  PT7,  PLA,  and  bleed  door  position  -  are  recorded  on  pencharts.  A  complete  list 
of  control  cab  parameters  is  also  given  in  (10). 

2.2  DUST  INJECTION  SYSTEM  (D5S). 

The  DIS  functions  to  deliver  a  continuous,  spatially  uniform  concentration 
of  dust  to  the  engine  face  at  velocities  representative  of  the  engine  flying  through  a 


quasi-stationary  dust  cloud.  It  is  designed  to  deliver  dust  concentrations  as  high  as 
1000  mg/m ^  having  sizes  up  to  250  microns. 

The  DIS  described  in  (10)  consists  of  a  compressed  air  supply,  a  weighbelt 
feeder,  a  feeder  controller,  a  dust  hopper,  a  set  of  venturis  and  nozzles,  and  an  engine 
inlet  tube.  Dust  falls  from  the  hopper  onto  the  belt  of  the  feeder.  The  speed  of  the 
belt  is  adjusted  by  the  controller  to  match  the  weight  flow  set  point  entered  in  the 
controller.  The  dust  then  falls  from  the  belt  and  is  carried  by  gravity  and  suction  into 
a  set  of  venturis  through  which  air  at  high  velocity  is  passing.  The  air  conveys  the 
dust  through  tubes  to  a  set  of  nozzles  at  the  intake  plane  of  the  engine  inlet  tube. 
The  dust  mixes  with  the  initially  clean  intake  air  as  it  travels  along  the  inlet  tube 
until  the  spatial  concentration  of  dust  is  uniform  at  the  engine  face.  Previous  use  of 
similar  hardware  indicated  a  set  of  six  nozzles  located  three  engine  inlet  diameters 
upstream  of  the  engine  injection  dust  at  the  velocity  of  the  clean  intake  air  would 
provide  sufficient  mixing.  Indeed  it  did;  the  intensity  of  the  glow  generated  by  the 
dust  striking  the  moving  fan  blades  varies  with  concentration  and  the  glow  intensity 
was  uniform.  The  exit  velocity  of  the  dust  is  set  by  metering  the  pressure  and  flow 
rate  of  the  conveying  air  upstream  of  the  venturi.  This  was  done  by  regulating  the 

compressed  air  supply  pressure  and  choking  the  flow  upstream  of  the  venturis. 

2.3  EXPERIMENTAL  PROCEDURE. 

All  engines  were  sent  through  a  set  of  preliminary  runs  for  the  purposes 
of  (a)  training  personnel  in  engine  operation;  (b)  checking  engine,  instrumentation,  and 
facility  operation;  and  (c)  generating  pre-exposure  engine  baseline  data.  Following  these 
runs,  the  engines  were  exposed  to  dust  at  concentrations,  power  settings,  and  exposure 
times  given  by  their  respective  test  matrices. 

The  test  matrices  were  based  on  knowledge  of  previous  engine  behavior 

and  environments  the  engines  would  be  expected  to  experience.  The  individual  runs 

were  chosen  (10-12)  to  give  an  indication  of  the  influence  of  concentration  and  power 
setting  on  damage  rate  and  damage  mechanism. 

Prior  to  each  exposure  run,  dust  weight  flow  was  calculated  to  match  the 
desired  dust  concentration.  The  volume  flow  of  air  was  inferred  through  the  use  of  a 
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nomograph  of  air  weight  flow  versus  power  setting  which  was  generated  during  the 
preliminary  runs. 

In  each  exposure  run,  the  engine  under  test  was  brought  to  the  EGT  set 
point  for  the  run,  allowed  to  thermally  stabilize,  and  then  exposed  to  dust  by  first 
turning  on  the  injecting  air  and  then  starting  the  weighbelt  feeder.  All  parameters 
were  recorded  just  before  dust  ingestion  and  several  times  each  hour  during  ingestion. 

Following  each  run,  the  last  stage  of  both  the  high  and  low  pressure 
compressors,  the  combustor,  and  the  front  and  rear  of  the  turbine  were  inspected  by 
borescope  for  signs  of  exposure  (e.g.  erosion,  glassification).  The  borescope  was  used 
most  frequently  to  get  a  qualitative  indication  of  the  status  of  the  tip  clearance  in 
the  high  pressure  compressor  because  early  in  the  testing  it  became  apparent  that  surge 
due  to  clearance  increases  would  be  the  dominant  engine  problem  for  the  class  of 
engines  being  used  here. 

2.4  ENGINE  DESCRIPTION. 

Results  obtained  from  two  engines,  the  P/W  TF33-P11A  and  the  P/W 
357-Pl*'V  are  described  in  this  report.  The  3 37  is  a  twin  spool  turbojet  engine.  At 
Takeof  Rated  Thrust  (TRT),  it  generates  approximately  12,000  lbs  of  thrust  while 
pumping  185  Ibs/sec  of  air  at  a  turbine  inlet  temperature  of  1640°F.  The  TF33  is  the 
turbofan  derivative  of  the  357.  It  has  one  less  low  pressure  compressor  stage  and  one 
more  iow  pressure  turbine  stage  than  the  357.  The  TF33  generates  18,000  lbs  of  thrust 
at  WRT,  pumps  460  lbs/sec  of  air,  and  has  a  turbine  inlet  temperature  of  1600°F.  It 
is  appropriate  to  run  these  engines  as  an  erosion  comparison  pair  because  the  two  have 
very  similar  internal  flow  geometries  and  their  HP  compressor  rotor  blades  are 
constructed  from  different  materials;  the  HP  compressor  blades  of  the  357-P19W  are 
stee'  and  those  of  the  TF33-P11A  are  titanium. 

2.5  OBSERVATIONS. 

Four  engines  have  been  run  in  the  LERC  in  dust  exposure  tests:  three 
TF33-Pl!As  and  one  357-P19W.  Dust  mixture,  time  at  assigned  power  settings  and 
concentration  have  varied  among  engines.  However,  after  several  hours  of  exposure, 
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the  response  of  all  engines  has  been  a  greater  susceptibility  to  surge  due  to  compressor 
erosion. 


After  nearly  eleven  cumulative  hours  of  dust  ingestion  testing  the  first 
engine  run  in  the  LERC,  a  TF33-P11A,  went  from  essentially  steady  state  operation  into 
a  series  of  surges  at  a  frequency  of  approximately  one  to  two  per  second.  The  condition 
was  relieved  by  retarding  the  throttle  to  idle.  During  the  hour  preceding  the  surges, 
the  engine's  exhaust  temperature,  pressure  ratio,  and  fuel  flow  steadily  dropped  with 
respect  to  high  pressure  compressor  speed.  Approximately  five  minutes  before  the 
surges  occurred,  the  engine  began  a  slight  autoacceleration;  fuel  flow,  pressure  ratio 
and  temperature  began  to  increase  without  any  operator  throttle  movement.  During 
the  surges,  the  pressure  ratio  rapidly  dropped  and  fluctuated  in  time  with  the  surges 
and  the  exhaust  gas  temperature  rose  by  over  300°F  in  ten  seconds. 

It  is  felt  that  this  sequence  of  events  occurred  because  of  the  inability 
of  the  eroded  high  pressure  compressor  to  sustain  the  control  set  point  pressure  in  the 
burner.  The  engine  was  operating  in  a  regime  in  which  compressor  discharge  pressure 
essentially  equals  burner  pressure.  The  trend  of  the  data  just  before  surge  suggests 
that  the  fuel  controller  increased  fuel  flow  to  bring  the  burner  to  its  set  pressure  and 
eventually  created  enough  backpressure  on  the  compressor  to  surge  it.  If  the  engine 
had  been  operating  at  a  higher  power  setting,  then  the  events  likely  would  not  have 
occurred  because  compressor  discharge  pressure  would  have  been  greater  than  burner 
pressure  and  surges  would  have  occurred  only  when  trying  to  accelerate  through  or 
operate  in  the  range  in  which  burner  and  compressor  discharge  pressures  are  essentially 
equal. 


Severe  compressor  erosion  accompanied  by  hard  engine  surges  describes 
the  behaviors  of  the  second  and  third  TF33  and  the  one  357.  Both  were  brought  into 
states  of  substantia!  deterioration  at  higher  power  settings  and  only  while  attempting 
to  accelerate  from  idle  did  surge  occur.  In  states  of  relatively  mild  deterioration, 
acceleration  rate  had  to  be  limited  and  after  heavy  deterioration,  the  engines  could  not 
be  accelerated  at  all.  It  also  became  apparent  with  these  engines  that  over  the  course 
of  exposure,  the  erosion  caused  a  readjustment  of  operating  parameters;  especially  rotor 
speeds  and  fuel  and  air  flows.  References  (10-12)  contain  a  complete  description  of 
the  observed  engine  behavior  as  a  function  of  dust  exposure. 
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SECTION  3 

EROSION  THEORY  FUNDAMENTALS 


Engine  components  which  are  exposed  to  dust  laden  air  sustain  damage 
by  continuous  pitting  and  cutting  of  the  metal  surfaces.  Impact  velocities  of  the 
particles  are  influenced  by  inlet  air  velocity  and  rotational  speed  of  the  rotor  blades. 
Typically,  the  impact  velocity  at  the  first  engine  rotor  is  600  ft/s  or  less.  After  the 
first  impact,  a  dust  particle— or  fragments  of  it— will  rebound,  may  accelerate  or 
decelerate  as  it  experiences  drag  in  the  air,  may  collide  with  other  particles,  and  then 
impact  again  on  moving  or  stationary  engine  surfaces. 

At  the  microscopic  level,  a  statistical  approach  must  be  used  in  describing 
the  erosion  process  because  the  simple  conceptual  case  of  a  single  particle  impacting  at 
a  given  velocity  and  angle  is  in  practice  "randomized"  by  several  complicating  factors. 
For  example,  the  particles  themselves  are  inhomogeneous  in  hardness  and  irregular  in 
shape  with  several  sharp  corners.  As  they  approach  the  metal  surface  with  no  preferential 
orientiation,  they  do  a  random  amount  of  work  on  it.  The  work  will  vary  from  very 
little— for  example  for  impact  broadside  on  a  flat  surface— to  transfer  of  a  substantial 
fraction  of  the  particle  momentum  into  deformation  of  target  material  or  shear  and 
impulse  of  target  fragments.  In  addition,  the  effective  angle  of  impact  is  randomized 
by  the  fact  that,  after  an  initial  incubation  period,  the  target  material  becomes  pitted 
with  craters  and  the  local  impact  angle  between  the  small  particle  and  the  eroded 
surface  may  deviate  considerably  from  the  approach  angle  before  the  impact. 
Furthermore,  the  particle  velocity  may  be  randomized  due  to  particle-particle  collisions, 
at  least  above  some  low  level  of  dust  concentration.  All  these  effects  have  been 
discussed  previously  in  the  literature  (13-17). 

A  methodology  to  predict  erosion  of  gas  turbine  components  must  separately 
define  the  particulate  environment  from  front  to  rear  of  the  engine,  the  shape  of  the 
target  area  and  the  associated  flowfield,  and  the  erosion  rate  for  the  special  target 
materials.  Although  in  an  engine  the  shape  of  the  target  surfaces  is  very  complex,  it  is 
in  fact  the  best  known— at  least  initially— among  the  above  factors.  The  flowfield  for 
the  uneroded  engine  is  also  complex,  but  known,  since  the  effect  of  the  presence  of 
the  particles  on  the  aerodynamics  of  the  engine  can  be  ignored  (to  a  first  approximation 
for  all  the  low-concentration  applications  of  interest  here).  On  the  other  hand,  local 
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changes  in  geometry  may  become  very  important  in  accelerating  or  decelerating  locally 
the  instantaneous  erosion  rate  as  the  engine  components  are  eroded. 


Among  the  fundamental  parameters  for  modeling  erosion,  the  erosion  rate 
is  conceptually  well  defined  and  may  be  obtained  experimentally  (18).  The  erosion  rate 
is  determined  by  measuring  the  weight  of  a  target  specimen  before  and  after  testing.  It 
is  expressed  as  milligrams  of  material  removed  per  gram  of  impinging  solid  particles. 
For  the  metallic  materials  of  interest  in  engines,  a  considerable  body  of  experimental 
data  exists.  However,  its  interpretation  is  not  always  clear  because  of  experimental 
differences  in  the  ways  the  data  were  obtained  and  the  still  inadequate  development 
of  erosion  theory.  For  example,  for  a  gas-particle  jet  impinging  on  a  target,  the 
amount  of  impinging  solid  particles  is  known  from  the  particulate  feed  rate  and  the 
exposure  time.  For  a  target  immersed  in  a  gas  particle  stream  a  capture  ratio  must 
be  determined.  This  is  the  case  for  the  moving  or  stationary  parts  of  an  engine  for 
which  the  cumulative  mass  of  dust  ingested  is  given  by: 

L  =  (particle  concentration)  x  (rate  of  volumetric  flow)  x  (exposure  time) 

where  all  quantities  are  taken  at  the  engine  inlet.  The  difficulty  in  determining 
appropriate  values  of  the  capture  ratio  at  specific  engine  locations  for  stator  vanes, 
rotor  blades  and  housing  surfaces  will  become  apparent  from  the  later  discussion  of 
gas-particle  dynamics  in  multistage  turbomachines. 

Limiting  considerations  to  well-defined  erosion  rate  data,  one  finds  that 
some  major  trends  are  well  established  and  bear  significantly  on  the  gas  turbine  erosion 
problem.  Experimental  observations  (16,  18,  19,  20)  indicate  that  dust  angle  of  attack 
and  velocity  strongly  influence  the  erosion  rate  for  a  given  combination  of  dust  and 
target  material.  Other  parameters  such  as  particle  size  and  shape,  and  dust  concentration 
may  have  a  lesser  effect,  if  any,  on  the  results. 

The  influence  of  angle  of  attack  on  the  removal  rate  of  material  is  one 
of  the  most  distinct  features  of  the  erosion  of  ductile  materials  such  as  c^g'n^  meta's. 
In  the  velocity  range  from  250  to  600  ft/s,  the  erosion  rate  increases  rapidly  with  angle 
of  attack  from  zero  at  0°  to  a  maximum  value  which  occurs  at  an  angle  of  approximately 
20°.  As  the  angle  of  attack  increases  further  from  this  value,  the  erosion  rate  decreases 
gradually  to  a  residual  value  at  90°.  A  relationship  of  the  erosion  rate  to  the  dust 

8 


velocity  based  on  the  supposition  that  the  erosion  process  be  proportional  to  the  kinetic 
energy  of  the  oncoming  particles  was  proposed  in  early  investigations.  However,  there 
is  a  lack  of  agreement  among  the  results  of  the  theories  based  on  kinetic  energy 
exchange  and  the  experimental  results.  This  is  still  a  perplexing  issue  in  the  field  of 
solid  particle  erosion.  The  conventional  expression  for  volume  removal  -  volume  removed 
per  particle  impact  -  is  given  as 

e  =  gvprybf  v  o) 

where  g  is  a  functional  representing  the  dependence  on  all  of  the  paramaters  not 
explicitly  stated  and, 

is  the  particle  impact  velocity 
rp  is  the  mean  radius  of  the  impacting  particle 

f  (  a  p)  is  the  function  of  the  dust  angle  of  approach  (discussed  above) 
a  is  the  velocity  exponent 
b  is  the  particle  size  exponent 

The  velocity  exponen*  has  been  generally  found  to  range  from  two  to 
three  for  materials  displaying  ductile  behavior  and  a  wide  range  of  particle  impact 
conditions;  however,  exponents  up  to  four  have  been  proposed  (6).  Note  that  an  even 
larger  variation  of  the  velocity  exponents  is  reported  in  the  literature  for  brittle 
materials  with  values  that  range  from  1  to  6.5  (21). 

Less  is  known  about  the  particle  size  exponent.  For  metals,  it  is  generally 
stated  that  the  volume  removal  is  independent  of  particle  size  for  particle  diameters 
in  excess  of  100  vm  (16).  However,  for  the  smaller  particle  sizes  which  are  of  interest 
in  turbomachinery  applications,  the  variations  in  particle  size  effects  may  be  significant. 
For  example,  in  brittle  materials  tested  with  particles  as  small  as  8  Pm,  the  particle 
size  exponent  was  found  to  vary  from  3  to  larger  than  4,  and  it  was  found  to  be 
coupled  with  the  velocity  exponent  as  well  (20). 

Erosion  in  muluotage  turbomachine^  is  significantly  more  complex  than 
the  erosion  of  simple  targets  such  as  flat  or  cylindrical  specimens.  This  is  due  primarily 
to  the  changes  which  the  erodent  particles  incur  from  the  inlet  to  the  exit  of  the 
engine  as  well  as  the  difficulties  in  tracing  the  particle  trajectories  through  the  engine 
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gas  path  after  multiple  impacts.  It  is  also  clear  that  the  cumulative  influence  of 
particle  dynamics  and  changes  in  the  dust  environment  makes  erosion  of  both  stationary 
and  moving  machine  components  very  dependent  on  their  stage  location  and  on  the 
specific  engine  design. 

The  motion  of  the  particles  which  are  convected  through  the  compressor 
and  the  turbine  should  be  analyzed  to  derive  bounds  on  the  location  and  severity  of 
component  erosion.  Considerations  related  to  the  trajectories  for  particles  having 
different  sizes  are  necessary;  (a)  to  estimate  the  erodent  weight  and  characteristics  at 
each  stage  location  and  (b)  to  estimate  heat-transfer  levels  while  the  particles  travel 
through  the  combustor  to  the  turbine  stages. 

Particles  suspended  in  a  gas  flow  move  under  the  influence  of  forces 
exerted  on  them  by  the  gas.  Ordinarily,  the  reaction  of  these  forces  affects  the  gas 
flow  only  to  a  negligible  extent  but  this  influence  may  have  to  be  considered  at  high 
concentration  or  for  highly  turbulent  flows.  More  importantly  for  the  engine  application, 
the  interaction  of  particles  with  wakes  of  other  particles  may  be  of  consequence.  The 
equation  of  motion  of  individual  particles  equate  the  product  of  mass  and  acceleration 
to  the  sum  of  various  forces;  namely,  viscous  drag,  gas  pressure  gradient  (or  alternatively 
gas  acceleration),  apparent  mass  acceleration,  and  Basset  force.  The  latter  force  accounts 
for  deviation  of  the  flow  pattern  around  the  particle  from  that  of  steady  flow.  In 
addition,  lateral  forces  on  the  particles  may  develop  as  a  result  of  particle  rotation, 
slip  in  the  shear  flow  of  a  boundary  layer,  and  other  effects.  Fortunately,  the  200:1 
ratio  between  the  density  of  the  particulate  and  the  density  of  the  gas  can  be  used 
to  simplify  the  analysis  in  the  application  of  interest  here.  With  such  a  high  density 
ratio,  it  can  be  shown  that  it  is  reasonably  justified  to  equate  the  mass-acceleration 
term  to  the  viscous  drag  only.  Because  the  flow  in  gas  turbines  is  severely  swirled, 
residual  effects  from  other  forces  cannot  be  ruled  out  entirely  and  the  conclusions 
derived  here  are  contingent  upon  experimental  verification. 


In  the  approximate  form,  the  equation  of  motion  for  a  spherical  particle 
convected  by  the  gas  becomes: 


(2) 
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where  the  substantial  derivative  D/Dt  is  taken  along  the  path  of  the  particle  and 

Dp  =  diameter  of  particle 

Up  =  velocity  of  particle 

u~  =  velocity  of  gas 

n,  pp  =  gas  and  dust  densities 

The  Reynolds  number  for  a  particle  can  be  expressed  as 

Rep  =  pDp  i  Cf  -  u^  I  (3) 

y 

In  our  application,  we  estimate 
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and  consequently 

Rep  =  0  to  2300 

For  this  range,  the  relationship  for  the  particle  drag  coefficient  can  be  given  as 

Cd  =  24  (1  +  1/6  Rep2/3)  (4) 


which  fits  the  standard  drag  curve  to  better  than  30%  and  the  trajectory  of  particles 
could  be  calculated  for  any  known  airflow  pattern  and  given  particle  size. 

A  useful  parameter  which  indicates  how  an  initial  velocity  difference  u 
-  Up  decays  in  time  in  a  constant  gas  velocity  is  easily  derived  from  Eq.  (3).  For  Rep 
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less  than  100  where  C<j  =  ^  to  a  good  approximation,  it  can  be  shown  that  an  initial 
velocity  difference  (u  -  Up)  §ecays  exponentially  with  a  relaxation  time  xv  given  by 

T  v  =  Dp  P  p 

18  u 

For  larger  Rep  values,  xv  remains  a  convenient  qualitative  measure  of 
velocity  relaxation.  For  the  present  application,  it  is  estimated  that  in  the  cooler 
(lower  viscosity)  gas  situation,  the  relaxation  times  are  in  the  range: 

(tv)l  =  0.03  ms  for  particles  1  um  in  diameter 

(  x  V)3Q  =  0.5  ms  for  particles  30  um  in  diameter 

The  above  relaxation  times  may  be  compared  with  the  time  required  for 
a  blade  to  transit  a  vane  exit  passage  in  order  to  determine  approximately  how  many 
times  a  particle  is  likely  to  strike  the  same  rotor  blade  and  where  the  impacts  are 
more  likely  to  occur.  However,  a  description  of  how  the  particles  will  rebound  upon 
impact  is  needed  first. 

Particle  rebound  correlations  were  obtained  by  Tabakoff  (22)  using  an 
experimental  approach  based  on  Laser  Doppler  Velocimeter  (LDV)  measurements.  A 
large  number  of  impacts  of  fly  ash  particles  having  a  diameter  of  15  um  were  observed 
as  they  impinged  on  metal  targets  at  320  ft/s.  Targets  of  annealed  2024  aluminum, 
6A1-4V  titanium  and  410  stainless  steel  were  used.  From  the  measurements,  values  3f 
the  average  velocity  restitution  ratio  (V2/V1)  and  the  average  directional  coefficient 
(  82/  8j)  were  obtained  for  impingement  angles  from  15  to  75°.  The  shape  of  the 
statistical  distributions  of  the  above  parameters  were  also  reported  because  they  are 
important  for  detailed  trajectory  calculations. 

The  velocity  restitution  ratio  parameter  is  a  measure  of  the  momentum 
lost  by  the  particle  at  impact  and  as  such  it  should  correlate  with  the  work  done  on 
the  target  material.  This  work,  in  turn,  determines  the  erosion  rate.  In  (22)  it  is 
shown  that  the  normal  velocity  restitition  ratio  Vjsj 2/^N  1  an^  the  tangential  restriction 
ratio  V-j-2/V-pj  do  not  contribute  equally  to  ductile  erosion.  In  fact,  the  average  value 
of  the  tangential  velocity  restitution  ratio  is  the  only  parameter  which  supports  the 

12 


theory  that  erosion  is  proportional  to  a  particle  momentum  loss  at  impact.  This  is 
because  the  maximum  erosion  of  the  metallic  materials  tested  occurs  at  values  of  the 
impingement  angle  where  the  tangential  velocity  restitution  ratio  is  low.  For  our 
purpose  the  generic  values  of  Table  1  were  used. 

Table  1.  Particle  rebound  characteristics. 
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It  should  be  noted  that  the  effect  of  approach  velocity  on  the  mean  value  of  the 
restitution  ratios  has  been  neglected  in  accordance  with  the  assumptions  and  findings 
of  several  authors  (23,  24). 

Using  the  above  particle  rebound  characteristics  and  the  relaxation  times 
computed  earlier,  the  frequency  and  location  of  impacts  of  a  particular  particle  on  a 
given  rotor  blade  were  estimated.  At  each  stage  it  was  assumed  that  all  trie  particles 
initially  impact  the  rotor  near  the  leading  edge.  The  results  of  this  very  simple  analysis 
indicate  that  nearly  all  of  the  particles  will  strike  the  same  rotor  Made  twice.  A 
small  fraction  of  the  particles  could  strike  the  blade  three  times.  The  blade  leading 
edge  and  trailing  edge  regions  experience  the  highest  number  of  impacts.  The  trend 
indicated  by  the  comparison  of  30  pm  particles  to  1  um  particles  is  that  erosion  of 
the  stationary  blades  will  decrease  and  erosion  of  the  rotating  blades  will  increase  if 
the  finer  particles  are  more  abundant.  In  general,  a  less  severe  erosion  is  indicated  in 
the  hot  turbine  flow  where  the  velocity  relaxation  times  are  considerably  smaller. 

In  (6)  and  (25),  detailed  predictions  of  particle  trajectories,  erosion  patterns 
and  eroded  mass  are  discussed.  The  results  were  obtained  using  a  sophisticated  Monte 
Carlo  simulation  model  of  the  erosion  process  in  relatively  simple  gas  turbine  cascades. 
A  one  and  one-half  stage  compressor  and  a  two  stage  turbine  case  were  considered. 
In  principle,  detailed  descriptions  of  erosion  intensity  and  pattern  can  be  obtained 
a  prior'  for  stationary  and  moving  gas  turbine  components  by  using  that  technique. 
However,  at  the  present  time  the  model  does  not  include  modeling  of  particle 
fragmentation  nor  does  it  account  for  erosion  rates  at  elevated  temperatures. 
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From  the  above  discussion  it  is  clear  that  a  number  of  fundamental  erosion 
issues  must  be  resolved  before  reliable  a  priori  predictions  of  erosion  in  realistic 
turbomachine  configurations  can  be  made.  At  present  the  two  most  important  areas 
for  which  basic  information  is  lacking  are  (a)  particle  breakup  models  and  d?a  and 
(b)  erosion  rates  at  elevated  temperatures. 

Having  recognized  this  situation,  the  remainder  of  this  study  followed  the 
approach  of  a  semiempirical  analysis  in  which  experimental  observations  relative  to  the 
level  and  pattern  of  erosion  in  the  TF33  and  357  engines  were  utilized  to  interpret 
the  operation  of  the  deteriorated  engines.  These  observations  could  also  be  used  to 
guide  in  the  prediction  of  deteriorated  engine  behavior  for  different  engine  types  which 
may  be  tested  in  the  future. 


SECTION  4 

FORMULATION  OF  PERFORMANCE  DETERIOR ATION 


Ingestion  of  dust-laden  air  has  the  effect  of  deteriorating  the  performance 
of  a  gas  turbine  engine  as  a  result  of  gas  path  component  erosion.  This  erosion  changes 
the  structural  character  and  thus  the  aeroelastic  behavior  of  the  engine  as  well  as  the 
aerodynamic  performance.  The  unsteady  aeroelastic  stability  is  not  expected  to  be 
rapidly  influenced  by  erosion.  However,  the  eroded  blades  are  known  to  quickly  lose 
the  ability  to  efficiently  increase  the  total  pressure  of  the  stage.  In  an  effort  to 
quantify  the  loss  in  performance  of  the  engine,  a  compressor  stage  will  be  discussed 
and  then  extrapolated  into  a  complete  compressor  map.  First  a  suitable  stage 
performance  map  will  be  calculated  from  the  pre-ingestion  compressor  map.  Then,  a 
deteriorated  stage  performance  will  be  determined  from  loss  correlations  and  measured 
blade  erosion.  The  deteriorated  stage  performance  maps  are  finally  stacked  to  arrive 
at  an  overall  deteriorated  compressor  map. 

To  estimate  the  change  in  stage  performance  with  respect  to  erosion,  the 
primary  sources  of  total  pressure  loss  affected  by  the  blade  aerodynamic  efficiency 
were  considered.  Attention  was  given  to  three  distinct  origins  of  total  pressure  loss: 
(a)  losses  associated  with  increasing  tip  clearance;  (b)  profile  losses  due  to  the  blade 
boundary  layer;  and  (c)  secondary  flow  losses  produced  at  the  hub  end  wall  and  casing 
boundary  layer  regions.  Of  considerable  interest  is  the  presence  of  total  pressure  losses 
in  the  tip  region  which  are  large  in  comparison  with  the  other  regions.  Within  the 
compressor,  the  impact  of  dust-laden  air  will  cause  considerably  greater  blade  leading 
edge  erosion  at  the  larger  span  locations  as  well  as  tip  region  erosion.  This  blade 
profile  modification  is  illustrated  in  Figure  I  which  illustrates  a  357  HP  compressor 
blade  that  has  been  subjected  to  a  dusty  environment  (11).  This  machine  uses  both  an 
LP  and  a  HP  compressor  and  the  erosion  of  both  components  is  described  in  (II).  It 
suffices  here  to  note  that  the  severity  of  erosion  increases  from  the  front  to  back 
stages  of  both  compressors,  being  most  severe  in  the  seven  stage  HP  compressor.  A 
region  of  significant  erosion  will  cause  a  total  pressure  distortion  in  the  radial  direction 
(26).  This  effect  leads  to  a  flow  instability  that  could  cause  boundary-layer  separation 
in  the  tip  region.  Such  behavior  will  degrade  the  stage  performance  and  may  lead  to 
a  compressor  surge. 
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4.1 


TIP  CLEARANCE  TOTAL  PRESSURE  LOSSES. 


Several  authors  (26-32)  have  presented  techniques  for  estimating  losses 
due  to  tip  clearance  leakage.  Experimental  results  (27)  show  a  linear  decrease  in 
efficiency,  with  increased  tip  clearance  to  blade  height  ratio,  Figure  2.  Theoretical 
models  have  shown  fair  agreement  (29)  with  experimental  results  in  certain  regions  of 
clearance  ratio. 


Recently  a  model  presented  in  (27)  was  found  to  agree  well  with  the 
experimental  result  of  seven  different  axial  flow  compressors  (see  Figure  2).  The  model 
accounts  for  leakage  loss  and  pressure  gradient  losses  associated  with  the  flow  interaction 
of  the  leakage  vortex.  The  suggested  expression  given  in  (27)  for  changes  in  efficiency 
in  an  axial  flow  compressor  blade  row  is  given  in  equation  5. 


A  n  _  .  7S9  s  *  j  q 

X  r/r.  c  sin '^3. 


1.17  <j> 


2  — 

a  r/r  tang 
t  m 
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7/rt  =  hub-tip  mean  radius  (r)/tip  radius  (rt) 
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The  variables  used  in  equation  5  are  local  values  at  the  blade  tip.  Figure  3  illustrates 
the  substantial  change  in  chord  and  thickness  experienced  in  the  tip  region  of  the  16th 
stage  HP  compressor  blade  of  the  357  machine  (II).  The  marked  erosion  at  the  leading 
edge  may  have  a  substantial  influence  on  the  pressure  distribution  on  the  blade  surface 
and  thus  result  in  a  decrease  in  the  strength  of  the  leakage  vortex.  For  this  reason, 
equation  4  is  expressed  in  terms  of  cascade  lift  coefficient  which  will  then  be  modeled 
for  changes  in  leading  edge  erosion.  The  definition  of  cascade  lift  coefficient  based  on 
mean  velocity  is  given  in  equation  6. 


c  =  ~  —  cos  e, 

[.  v  i  m 


(6) 


where  <P  is  pressure  coefficient 
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Figure  2.  Decrement  of  efficiency  due  to  tip  clearance  ratio  seven  axial  impellers. 
(Taken  from  [1] .) 


s/c 


is  flow  coefficient 

is  solidity  at  the  midspan 

is  the  mean  fluid  angle  at  midspan 


3 


m 


5tiP 


.759 
r/  r 


a  C, 


1.17  <J> 


3 

sm  6  cos  0 
m  m 


2  — 

a  r/r  tan  3 

t  u  J 


ip  X 
is 


(7) 


The  definition  of  efficiency  is  substituted  to  obtain  an  expression  for  the  total  pressure 
loss  coefficient  due  to  tip  clearance  effects.  Equation  7  can  be  used  to  estimate  the 
total  pressure  loss  induced  by  increased  tip  clearance  and  by  changes  in  the  lift 
coefficient  resulting  from  blade  profile  changes  at  the  tip. 

As  illustrated  earlier,  the  principal  effects  of  erosion  on  the  airfoil  are 
removal  of  material  in  the  pressure  surface  leading  edge  region,  thinning  of  the  trailing 
edge,  and  reduced  leading  edge  radius.  Reference  7  shows  a  significant  increase  in 
two  dimensional  profile  losses  as  a  result  of  blade  erosion.  Due  to  the  high  loading  at 
the  leading  edge  of  an  airfoil,  it  is  expected  that  any  change  in  this  region  will  be 
significant.  The  effect  of  thickness  distribution  along  the  airfoil  was  assumed  to  have 
a  minor  effect  on  the  overall  aerodynamic  performance  of  the  airfoil.  The  effect  of 
leading  edge  erosion  will  be  modeled  by  a  change  in  the  mean  camber  line  of  the 
airfoil  (see  Figure  4).  This  change  will  result  in  an  increase  in  the  angle  of  zero  lift 
and  thus  a  decrease  in  blade  lift  coefficient  at  fixed  inlet  flow  angle.  A  reasonable 
estimate  of  the  contribution  to  the  lift  coefficient  from  the  mean  camber  line  is 
accomplished  using  thin  airfoil  theory.  It  is  demonstrated  in  Appendix  A  that  for  any 
given  eroded  state, 

A  C?  =1.5  *  E  (  AC  )  (8) 

Equation  8  approximates  the  change  in  two  dimensional  blade  lift  coefficient  due  to 
leading  edge  erosion  and  is  dependent  on  the  height  of  the  uneroded  camber  line  at 
the  eroded  leading  edge.  Note  that  the  change  in  blade  lift  coefficient  from  equation 
8  should  be  subtracted  from  the  uneroded  blade  lift  coefficient  when  used  to  estimate 
profile  losses.  For  the  purpose  of  this  report,  the  uneroded  blade  lift  coefficient  is 
obtained  by  a  method  outlined  in  the  following  paragraph. 


Figure  3.  Blade  profile  erosion  near  the  tip. 


Figure  4.  Variation  of  mean  camber  line  with  erosion. 


It  is  not  recommended  that  classical  thin  airfoil  theory  be  used  to  find 
the  uneroded  blade  lift  coefficient  since  turbomachinery  blades  are  typically  highly 
cambered.  Rather,  an  improved  potential  flow  solution  was  obtained  by  using  a 
modification  of  classical  airfoil  theory  which  removes  the  assumptions  of  small  angle 
of  attack  and  small  camber  angle.  The  method  used  in  this  development  is  given  in 
Appendix  B  and  does  not  require  extensive  computational  facilities.  The  result  for  the 
357  profile  is  given  below. 

=  4.724  sin  a  +  1.614  cos  a  (9) 

4.2  MODELING  BLADE  PROFILE  LOSSES. 

Thus  far  the  development  has  concentrated  on  the  tip  region.  The 
considerations  made  in  deriving  the  local  apply  also  to  any  other  spanwise  location. 
In  order  to  relate  the  blade  lift  coefficient  to  the  cascade  lift  coefficient,  a  modifying 
coefficient  (described  in  (34)  and  (35)),  K,  must  be  used  to  account  for  the  effects  of 
channel  flow  in  the  cascade.  Although  the  use  of  a  cascade  coefficient  may  seem  to 
oversimplify  the  formulation  of  the  cascade  lift  coefficient,  it  has  been  found  to  be 
reasonable  in  subsonic  flow  (34,35).  Typical  values  of  K  are  approximately  0.5  for 
axial  flow  compressors.  For  small  changes  in  chord,  K  can  be  assumed  to  remain 
constant. 


Now  that  a  means  of  tracing  changes  in  cascade  lift  coefficient  with 
respect  to  leading  edge  erosion  exists,  a  method  of  relating  lift  coefficient  to  cascade 
profile  losses  must  be  developed.  From  a  two-dimensional  control  volume  analysis  on 
a  cascade,  expressions  for  lift  and  drag  coefficients  have  been  formulated  (36)  and  are 
given  below. 


Cd  -  (AP^l  s  cos  g  m  (10) 
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(11) 


r  _  9  --  cos  gm  (tan  Sj  -  tan  B  iA)  -  Cp  tan  Bm 
_  z  c 


The  loss  coefficient  based  on  mean  velocity  is  given  by  the  total  pressure  loss  across 
the  blade  divided  by  the  dynamic  pressure. 


PA 

1/2  „  w* 


(12) 


Equation  1 1  was  developed  from  a  two  dimensional  analysis  and  may  be  solved  for  drag 
coefficient  to  arrive  at  a  two-dimensional  profile  loss  coefficient  with  the  use  of 
equations  11  and  12.  To  express  loss  coefficient  with  respect  to  tip  rotational  speed, 
an  assumption  of  constant  axial  velocity  through  the  blade  row  is  made  to  yield 
equation  13. 
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There  are  numerous  detailed  accounts  of  two-dimensional  cascade  tests  which  will  give 
a  good  estimate  of  the  profile  losses  (37).  With  the  stage  performance  map  given  and 
the  blade  lift  coefficient  computed,  where  a  =  gj  -  v  ,  the  modifying  coefficient  may 
be  obtained  using  equation  5. 

In  addition,  from  knowledge  of  uneroded  lift  coefficient,  profile  loss 
coefficient,  and  blade  row  geometry  the  tan  gj  -  tan  g  jj  term  may  be  computed. 
The  blade  profile  of  a  J57  16th  stage  compressor  blade  was  substituted  into  equations 
A-9  and  A-4  to  yield  changes  in  lift  and  total  pressure  loss  coefficient  with  respect  to 
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leading  edge  erosion,  (see  Figures  5  and  6).  Since  erosion  is  concentrated  in  the  tip 
region,  the  integrated  profile  loss  increase  is  not  represented  by  Figure  6.  Equation 
15  accounts  for  the  uneroded  portion  of  the  blade  near  the  hub  and  provides  the 
averaged  total  pressure  loss  coefficient  due  to  leading  edge  erosion  in  the  tip  region. 
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where 


ue  =  fraction  of  blade  height  influenced 
4.3  SECONDARY  FLOW  LOSSES. 

The  losses  due  to  secondary  flows  must  also  be  examined  due  to  their 
strong  dependence  on  lift  coefficient.  The  literature  contains  correlations  (e.g.  34)  of 
losses  due  to  secondary  flows  given  by  Cq  =  .018  or  referring  to  equations  9  and  13 
a  loss  coefficient  given  by  equation  14. 


-sec  =  -018  <i>2  c/s  cos  Bm  (16) 

Since  the  secondary  flow  is  the  result  of  the  hub  endwall  and  casing  boundary  layers 
and  their  respective  local  lift  coefficients,  equation  16  will  be  separated  into  two  parts. 
The  correlation  will  be  assumed  to  have  equal  contributions  from  endwall  and  casing 
boundary  layers  for  the  uneroded  condition.  As  the  tip  region  erodes  the  local  two 
dimensional  lift  coefficient  will  decrease  and  thus  secondary  flow  losses  in  the  tip 
region  decrease.  The  contribution  due  to  the  hub  endwall  region  will  be  assumed 
constant  due  to  relatively  little  erosion  in  this  region.  The  new  expression  for  secondary 
flow  losses  given  in  equation  17  and  may  be  calculated  from  the  previous  formulation 
of  lift  coefficient. 
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4.4  STAGE  PERFORMANCE  DETERIORATION. 

In  order  to  study  the  effect  of  erosion  on  overall  compressor  deterioration, 
a  loss  in  stage  performance  will  be  modeled.  Experimental  studies  presented  by  (26) 
show  the  effect  of  tip  leakage  flow  on  a  stage  performance  map  (Figure  7).  The  study 
indicates  that  a  change  in  slope  does  not  occur,  but  rather  an  overall  shift  of  the 
pressure  coefficient  and  efficiency  versus  flow  coefficient  curves  produces  the  observed 
result.  This  behavior  will  allow  the  changes  in  performance  to  be  modeled  at  +he 
design  point  and  then  used  to  produce  an  entire  eroded  performance  map.  This 
assumption  is  of  particular  importance  since  the  modeling  results  are  not  expected  to 
yield  accurate  predictions  at  off  design  operating  points.  To  estimate  efficiency  and 
performance,  the  previous  total  pressure  loss  formulations  will  be  used  with  the 
knowledge  of  the  stage  characteristics.  Therefore,  stage  to  stage  fluid  angle  must  be 
known.  If  the  stage  flow  characteristics  are  not  available,  as  is  the  case  for  many 
older  engines,  an  estimate  of  fluid  angles  must  be  made  from  a  given  stage  performance 
map.  Compressor  stages  are  typically  designed  to  have  approximately  a  fifty  percent 
stage  reaction  for  peak  efficiency.  The  flow  coefficient  may  be  related  to  the  relative 
mean  fluid  angle  by  equation  18. 

4>  =  1/2  cot  gm  (18) 


The  design  point  is  taken  to  be  at  a  constant  flow  coefficient  as  the  stage  deteriorates. 
On  this  basis,  the  mean  fluid  angle  and  difference  in  tangents  of  inlet  and  exit  fluid 
angles  will  be  assumed  constant  at  the  design  point.  The  assumption  of  fifty  percent 
stage  reaction  also  gives  an  equal  amount  of  losses  through  rotors  and  stators.  The 
computed  and  measured  stage  characteristics  are  given  in  Table  2  for  the  13th  stage 
of  the  357  HP  compressor. 
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Table  2.  13th  stage  characteristics. 


mean  to  tip  radius 
solidity 

design  flow  coefficient 
design  pressure  coefficient 
tip  clearance 
blade  height 
mean  fluid  angle 

stagger  angle 
polytropic  efficiency 
profile  loss  coefficient 
cascade  modifying  coefficient 
cascade  lift  coefficient 
total  loss  coefficient 


r/rt 

0.92 

c/s 

1.0 

* 

0.685 

P 

0.543 

T 

0.05” 

h 

2" 

gm 

36° 

tan  8i  -  tan  Bj; 

0.464 

V 

30° 

n  pc 

0.86 

£  profile 

0.012 

K 

0.20 

cL 

0.64 

K 

0.050 

To  estimate  the  inlet  relative  fluid  angle,  and  thus  angle  of  attack  to  be  used  in 
equation  18,  a  typical  velocity  triangle  at  the  inlet  was  constructed  (Figure  8).  Equation 
19  was  used  to  estimate  the  inlet  fluid  angle. 

8  j  =  tan-1  (1/  4>)  -  R  (19) 

where  e  *  5 


The  measured  blade  erosion  and  calculated  changes  in  pressure  roefflei^nt  and  efficiency 
at  the  design  point  are  given  in  Table  3  for  the  13th  stage. 


Table  3.  Performance  change  due  to  erosion,  13th  stage. 


chord  change  at  tip 

A^tc 

0.06 

average  chord  change 

Ac/c 

0.0.3 

leading  edge  influence 

we 

0.3 

decrease  in  efficiency 

Anst 

0.064 

decrease  in  pressure  coefficient 

A\p 

0.041 

tip  clearance  increase 

At 

0.05" 

increase  in  loss  coefficient 

1.39 

The  above  calculations  were  made  assuming  no  increases  in  losses  through 
the  stator.  The  distribution  of  total  pressure  losses  for  the  three  contributors  discussed 
earlier  is  given  in  Table  4  below.  Formulation  of  loss  in  performance  was  carried  out 
for  all  seven  stages  of  the  high  pressure  compressor  and  are  presented  in  Table  5. 
These  loss  estimates  were  computed  with  the  use  of  the  stage  performance  formulated 
from  a  69.9  percent  speed  line,  which  will  be  discussed  later. 
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Table  4.  Distribution  of  total  pressure  losses  in  the  13th  stage. 


CLEAN 

O, 

0 

ERODED 

0. 

0 

Profi  le 

0.0125 

1 

' 

24 

0.0212 

24 

Secondary  flows 

0.0035 

6 

0.0027 

4 

Tip  clearance 

0.0340 

70 

0.0667 

74 

Total 

0.050 

100 

0.091 

100 

Table  5.  Computed  performance  loss  in  stages. 


STAGE 

Aiji 

Anst 

Z/K  total 

0 

10 

— 

0.010 

0.018 

1.10 

u 

0.036 

0.062 

1.47 

12 

0.047 

0.079 

1  .51 

13 

0.050 

0.084 

1.51 

14 

0.113 

0.  190 

2.13 

15 

0.120 

0.200 

2.14 

16 

0.139 

0.232 

2.23 
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SECTION  5 

MODELING  COMPRESSOR  RESPONSE 

5.1  FORMULATING  AN  AVERAGED  STAGE  PERFORMANCE  MAP. 

The  procedure  described  in  the  previous  section  can  be  used  to  determine 
stage  performance  deterioration  when  the  actual  stage  performance  characteristics  are 
given.  This  was  not  the  case  for  the  TF33  and  3  57  engines  which  were  used  in  these 
measurements.  Therefore,  a  method  was  developed  to  obtain  the  stage  performance 
map  from  the  overall  compressor  map.  It  is  important  to  incorporate  a  compression 
efficiency  in  the  method  so  that  off-design  performance  can  be  computed.  The  only 
available  input  parameters  are  assumed  to  be  the  compressor  map  with  corresponding 
overall  adiabatic  efficiencies.  The  method  outlined  below  is  based  on  the  relationship 
between  a  constant  rotational  speed  line  and  the  stage  total  pressure  coefficient  curve. 
The  compressor  map  provides  the  relationship  among  mass  flow,  total  pressure  rise, 
and  overall  efficiencies  along  constant  rotational  speed  lines.  The  inlet  total  pressure 
was  assumed  to  be  atmospheric  pressure  and  the  density  was  choosen  to  be  that  at 
standard  conditions.  The  exit  velocity  was  calculated  knowing  the  inlet  conditions, 
compressor  flow  areas,  and  the  efficiency.  The  total  pressure  increase  can  be  averaged 
by  the  number  of  stages  and  a  pressure  coefficient  for  inlet  and  exit  stages  can  be 
averaged.  A  suitable  averaged  stage  map  can  then  be  used  in  a  stage  stack-up  procedure. 
Each  stage  was  assumed  to  have  the  same  performance  map  but  different  operating 
points.  The  exit  conditions  can  then  be  calculated  from  the  continuity  equation  and 
the  definition  of  total  pressure,  equations  20  and  21. 


=  p4  V4 


(20) 


l/: 


(21) 
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In  order  to  obtain  a  third  equation  and  thus  allow  the  calculation  of  density,  velocity, 
and  pressure  at  the  exit,  the  second  law  of  thermodynamics  'and  the  definition  of 
polytropic  compression  efficiency  were  used  (see  Appendix  C). 

The  average  flow  coefficient  is  computed  from  the  average  of  inlet  and 
exit  velocity  divided  by  tip  rotational  speed.  The  average  pressure  coefficient  is 
computed  from  the  averaged  total  pressure  rise  per  stage  divided  by  the  reduced  density, 
equation  22. 
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By  using  the  reduced  density,  the  method  essentially  averages  the  input  mechanical 
energy  per  stage.  Figure  9  shows  the  computed  averaged  stage  performance  map  for 
a  69.9  percent  rotational  speed  line.  This  map  was  used  to  reconstruct  the  constant 
speed  line  of  the  compressor  map  by  e  stage  stacking  technique  described  later.  A 
slight  adjustment  was  needed  to  reproduce  the  compressor  map  line.  This  was  due  to 
the  lack  of  available  compressor  efficiency  data  which  required  that  estimates  of  local 
values  be  made.  From  the  previous  calculations  of  changes  in  stage  performance  an 
eroded  performance  map  for  several  stages  is  given  in  Figure  10.  If  stage  performance 
data  were  available,  they  would  be  preferable  to  an  averaged  stage  map.  Also,  the 
above  averaged  performance  map  will  only  yield  accurate  results  for  the  constant  speed 
line  used. 


5.2  GENERATING  AN  ERODED  COMPRESSOR  MAP  FROM  STAGE 

STACKING. 

The  compressor  map  foi  an  eroded  compressor  can  be  obtainec  ^>y  reversing 
the  procedure  described  in  Section  >.l  and  by  using  blade  erosion  data.  In  other  words 
the  deteriorated  stages  are  "stacked"  to  obtain  the  eroded  compressor  map.  Figure  1 1 
shows  the  clean  69.9  percent  constant  speed  line  and  a  corresponding  eroded  line 
computed  using  this  stage  stacking  technique.  It  is  assumed  that  the  operating  line  of 
the  compressor  does  not  shift  due  to  the  relatively  unchanged  fuel  controller,  combustor, 
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POLYTROPIC  EFFICIENCY  VS.  FLOW  COEFFICIENT 


Figure  9.  Estimated  stage  performance  map  from  the  69.9  percent  N2  line. 


POLYTROPIC  EFFICIENCY  VS.  FLOW  COEFFICIENT 


/APt  \ 

\v,pujj 


PRESSURE  COEFFICIENT  VS.  FLOW  COEFFICIENT 


0  = 


U* 


Figure  10.  Estimated  stage  performance  deterioration  from  measured  erosio 
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and  turbine.  The  shift  in  speed  lines  does  illustrate  a  decrease  in  corrected  weight 
flow  and  total  pressure  rise  at  a  constant  rotational  speed.  The  loss  in  corrected 
weight  flow  at  the  constant  speed  line  of  69.9%  was  computed  to  be  4.4%  of  the 
expected  weight  flow  of  the  undeteriorated  engine  at  the  operating  line. 

To  examine  the  change  in  surge  line  with  respect  to  erosion,  a  comparison 
of  compressor  stability  will  be  made  between  an  undeteriorated  and  a  deteriorated  low 
speed  line.  For  constant  properties,  the  inlet  axial  velocity  is  given  by  the  weight 
flow  through  the  compressor.  Therefore,  the  inlet  flow  coefficient  will  be  the  same 
for  both  the  deteriorated  and  undeteriorated  speed  line  at  the  weight  flow  corresponding 
to  the  undeteriorated  compressor  surge  line.  In  the  case  of  a  low  speed  surge,  the 
front  stages  are  at  high  incidence,  due  to  lower  flow  coefficient,  and  the  back  stages 
are  approaching  a  choking  condition.  With  increasing  stage  deterioration,  the  back 
stages  (at  the  same  speed)  are  experiencing  a  lower  density  relative  to  their  uneroded 
operation.  This  will  cause  the  exit  velocity  to  be  elevated  even  higher  than  the 
uneroded  operating  value.  The  deteriorated  engine  then  moves  closer  to  choking  the 
rear  stages  with  the  same  inlet  flow  coefficient.  With  the  flow  coefficient  corresponding 
to  stage  stall  unchanged  for  the  deteriorated  stage,  the  front  stages  remain  at  their 
unstable  high  incidence  values.  It  is  reasonable  to  assume  then  that  compressor  surge 
will  occur  at  the  same  corrected  weight  flow  as  that  of  the  undeteriorated  compressor 
at  low  speed  operation.  The  loss  in  available  surge  margin  and  surge  pressure  ratio  at 
the  69.9%  speed  surge  point  were  estimated  by  equations  23  and  24  and  are  presented 
in  Table  6  (38).  (See  Figure  12).  The  reduction  in  surge  margin  is  predicted  to  be 
51%  while  the  loss  in  surge  pressure  ratio  is  estimated  to  be  45%. 


ASM  =  PRS  -  PR0 
PR^ 


(23) 


Aprs  =  prs  -  pr0 

~pr7~' 


(24) 
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Figure  12.  Illustration  of  surge  margin  computation. 
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Table  6.  Loss  in  compressor  stability. 


Uneroded 

Eroded 

%  Decrease 

ASM 

.305 

00 

•a- 

• 

51% 

APRS 

.234 

.129 

45% 

5.3  EFFECT  OF  LPC  AND  HPC  COUPLING. 


In  matching  the  components  of  a  two-spool  gas  turbine,  it  is  convenient 
to  match  the  components  of  the  HP  spool  and  then  match  the  HP  spool  with  the  LP 
spool  components.  In  this  discussion,  the  influence  of  the  high  pressure  compressor 
(HPC)  deterioration  on  the  low  pressure  compressor  (LPC)  is  described.  The  matching 
map  for  the  HP  spool,  which  controls  the  HP  compressor,  is  altered  while  sustaining 
the  same  operating  line.  A  direct  method  of  matching  LP  and  HP  spools  is  described 
in  (39)  and  is  outlined  below  in  order  to  illustrate  the  qualitative  aspects  of  the  influence 
of  deterioration  on  spool  matching.  Due  to  the  aerodynamic  coupling,  only  two 
independent  conditions  are  needed  to  determine  the  operating  points  of  the  components. 
An  operating  point  is  assigned  to  the  HP  spool  and  the  LP  spool  matching  is  achieved 
directly.  A  schematic  of  the  numbering  system  used  is  given  in  Figure  13.  The  power 
and  continuity  relations  are  expressed  with  regard  to  the  LPC  and  LPT,  equations  25 
and  26. 
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Figures  14  and  15  represent  LPC  and  LPT  matching  maps,  respectively.  To  match,  an 
inner  spool  operating  point  is  taken  and  values  ol  W7/g"?/62,W3/93/6[.  ,  and  P^/P., 
are  noted  on  Figure  16.  Figures  14  and  15  can  be  superposed  so  as  to  satisfy  equations 
25  and  26.  The  ordinates  in  Figures  14  and  15  are  offset  by  the  value  of  the  accessory 
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Figure  15.  Low  pressure  turbine  matching  map. 
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Figure  16.  High  pressure  compressor  matching  map. 
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power  term  in  equation  25  and  the  abscissas  are  offset  by  the  calculated  value  of  the  log 

term.  The  operating  points  of  the  LPC  and  LPT  corresponding  to 
the  selected  HP  spool  operating  point  are  at  the  intersection  of  the  appropriate  lines 
of  constant  W?/oo/52  and  W5/e5/S5. 

To  evaluate  the  influence  of  HPC  deterioration,  the  lines  of  constant  N2 
are  shifted  down  in  Figure  16  and  the  above  procedure  is  reDeated.  The  LPC  and  LPT 
are  assumed  unchanged  so  their  matching  maps  are  not  altered.  Note  the  operating 
line  of  Figure  2  is  assumed  constant  so  a  decrease  in  W2*/,02//^2’^S>/,®5^B  anc*  ^5^2 
is  present  at  a  constant  inner  speed.  The  accessory  power,  fuel  air  ratio,  and  bypass 
ratio  terms  are  assumed  constant.  The  decreasing  P5/P2  value  will  decrease  the 
difference  between  the  abscissas  of  Figures  14  and  15.  The  ordinates  are  unchanged 
due  to  constant  accessory  power.  Figure  17  shows  the  comparison  between  the 

undeteriorated  and  deteriorated  by  superposition  of  Figures  14  and  15.  For  a  constant 
HP  spool  speed,  the  result  is  a  decrease  in  the  LP,  spool  speed  for  deteriorated 
operation.  Since  the  LPT  and  HPC  operating  line  is  assumed  unchanged  the  operating 
line  of  the  LPC  may  be  assumed  constant.  This  is  a  valid  assumption  for  either  a 
small  shift  in  the  operating  line  or  an  operating  point  which  is  moving  along  an  unchanged 
operating  line.  The  change  in  weight  flow  through  the  LPC  is  given  by  the  continuity 
relation  of  equations  27  and  28. 

AW3 

AW1  =  (1-B9) (I-B3)  (27) 


AW1 


(28) 


The  deteriorated  LP  spool  speed  may  be  found  with  the  appropriate  figures  for  the 
given  engine. 
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Hgure  17.  Superposition  of  LPC  and  LPT  matching  maps;  deteriorated  and  undeteriorated. 
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SECTION  6 

DETERIORATED  ENGINE  RESPONSE 

6.1  BLEED  AIR  CONTROLS. 

Given  the  357-P19W  and  TF33-P11A  engines,  there  are  two  ways  the 
engine  operator  can  alter  the  compressor  inlet  and  exit  conditions  to  make  them  less 
prone  to  surge.  They  are  the  intercompressor  bleed  and  the  anti-icing  air  systems. 

The  intercompressor  bleed  is  essentially  a  valve  actuated  according  to  a 
pressure  schedule  which  allows  a  portion  of  the  air  between  the  two  spools  to  dump 
overboard.  The  valve  remains  open  at  low  speeds,  allowing  the  LPC  to  essentially 
freewheel  (i.e.,  accelerate  gas  but  compress  it  very  little).  As  the  engine  is  accelerated 
through  a  scheduled  LPC  pressure  ratio,  the  valve  closes  and  the  spools  rematch. 
Pressure  ratio  across  the  LPC  rapidly  rises  and  the  HPC  pressure  ratio  first  falls  and 
then  rises  with  increasing  thrust. 

By  being  presented  with  air  at  lower  density  and  velocity  than  if  the  valve 
were  closed  with  air  at  lower  density  and  velocity  than  if  the  valve  were  closed,  the 
HPC  is  less  likely  to  experience  the  low-speed  front-stage  surge  generated  by  high 
angle  of  attack  from  rear  stage  choking.  When  the  valve  does  close  in  normal  operation, 
density  and  velocity  at  the  high  compressor  rapidly  increase  but  at  this  rotor  speed 
the  compressor  cannot  sufficiently  compress  the  air  so  as  to  not  choke  the  rear  stages. 
If  the  HP  spool  had  been  turning  slower  than  expected  at  closing,  it  would  have  been 
more  likely  to  surge. 

The  anti-ice  air  system  serves  to  raise  the  temperature  of  the  metal 
surface  at  the  face  of  the  engine  in  order  to  avoid  potentially  dangerous  ice  buildup. 
The  system  routes  hot  air  from  the  cowling  downstream  of  the  HPC  through  a  thermally 
adjusted  valve,  inside  the  first  fan  stator  and  to  the  engine  inlet.  The  valve  limits 
the  anti-ice  flow  to  only  the  amount  necessary  to  be  effective  against  the  ambient 
conditions.  The  amount  bled  is  at  most  four  or  five  percent  of  the  total  flow.  By 
bleeding  air  from  the  exit  end  of  the  HPC,  the  anti-ice  air  system  effectively  reduces 
the  flow  restriction  felt  by  the  front  stages  and  moves  it  away  from  surge  at  all 
speeds;  particularly,  however,  the  front  stages  at  low  speed. 
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Both  of  the  systems  described  above  were  used  during  the  measurement 
program  to  make  acceleration  to  higher  power  settings  possible.  They  were  found  to 
become  effective  at  different  deterioration  and  power  levels  than  each  other.  One  of 
the  TF33s  eroded  to  a  point  at  which  it  could  not  be  accelerated  through  a  power 
setting  just  above  idle  unless  the  anti-ice  air  system  was  operating.  The  power  setting 
involved  was  low  enough  for  the  bleed  door  to  be  in  the  open  position  and  any  surge  in 
this  regime  would  be  in  the  high  pressure  compressor  front  stages.  As  the  compressors 
eroded,  the  largest  proportionate  material  loss  was  in  the  rear  stages.  As  described 
in  Section  5,  severe  erosion  in  the  rear  stages  results  in  lower  densities  and  higher 
velocities  for  a  given  speed  and  hence  even  more  of  a  tendency  to  choke  at  low 
rotational  speed. 

The  amount  of  air  bled  for  anti-icing  is  scheduled  based  on  ambient  inlet 
air  temperature  with  larger  amounts  bled  at  lower  temperatures.  This  method  of  surge 
avoidance,  then,  may  not  be  effective  at  high  ambient  temperatures  unless  the  thermally 
controlled  valve  were  overridden  (which  is  not  a  feature  of  these  engines).  Anti-ice 
air  did  not  have  a  noticeable  influence  on  surge  margin  other  than  at  the  low  power 
settings  previously  mentioned.  Rapid  accelerations  from  midpower  to  high  power  levels 
with  anti-ice  air  both  on  and  off  gave  essentially  identical  results. 

After  further  exposure,  the  TF33  experienced  surges  at  steady  state  just 
above  the  bleed  door  closing  point  and  a  single  surge  each  time  it  was  accelerated 
through  closing.  The  HPC  is  operating  at  nearly  its  design  optimum  speed  at  the  time 
of  bleed  door  closing  and  is  relatively  difficult  to  send  into  surge.  The  LPC  is  at  the 
same  time  rotating  at  low  speed  with  relatively  little  margin  and  is  likely  to  go  to 
front  stage  surge. 

As  erosion  progressed,  the  engine  spools  rematched  to  reflect  the  change 
in  work  distribution  between  them.  Typically,  the  air  weight  flow  for  a  given  engine 
speed  and  the  low  compressor  speed  for  a  given  EPR  decreased.  High  compressor 
speed  for  a  given  EPR  increased.  That  both  low  compressor  speed  and  air  flow  decreased 
is  to  be  expected  in  an  undamaged  compressor  and  indeed  the  low  spools  from  the 
engines  were  relatively  undamaged.  The  reason  the  high  compressor  pumps  less  at 
higher  speed  is  examined  in  Section  5.  The  role  of  the  fuel  controller  in  the  rematch 
is  not  clear  at  this  time. 


The  fuel  controller  sets  fuel  flow  based  on  a  preset  schedule  of  HPC 
speed,  burner  pressure,  and  power  lever  angle  (PLA).  A  sketch  of  a  fuel  controller  map 
is  shown  in  Figure  18.  For  a  given  PLA,  the  controller  maintains  a  fixed  fuel  flow- 
to-burner  pressure  ratio  and  a  fixed  compressor  speed  by  metering  fuel  to  the  combustor. 
Whenever  the  power  setting  is  above  cruise,  burner  pressure  is  constant  and  the  controller 
adjusts  fuel  flow  to  maintain  control  speed.  For  power  settings  below  cruise,  fuel 
flow,  burner  pressure,  and  speed  all  increase  with  increasing  power.  Fuel  flow  increases 
faster  than  burner  pressure  with  increasing  speed. 

It  is  felt  the  fuel  controller  contributed  to  the  quasi-steady  state  surge 
described  in  Section  2.7.  During  the  hour  preceding  surge,  the  PLA  was  increased 
several  times  to  compensate  for  decreasing  EPR.  HPC  speed  had  remained  constant 
and  both  fuel  flow  and  burner  pressure  had  been  falling.  Each  time  the  PLA  was 
adjusted,  burner  pressure  increased  to  its  "pre-fair’  value  and  fuel  flow  exceeded  its 
value.  On  the  fuel  control  map  this  is  represented  by  a  leftward  shift  of  the  constant 
PLA  lines  and  a  new  operating  point  above  the  operating  line.  The  fuel  controller 
response  can  be  either  to  cut  fuel  flow  (which  it  did  several  times  initially)  or  to 
search  for  a  new  one  line  operating  point  at  a  higher  speed  (which  it  did  before  going 
into  surge).  This  would  mean  the  instantaneous  operating  point  would  move  to  the  left 
and  the  fuel  flow  would  begin  to  increase  to  return  to  the  line.  If  this  increased 
burner  pressure  sufficiently  to  overpower  the  exit  pressure  of  the  HPC,  then  surge 
would  occur. 

6.2  DETERIORATION  IN  WEIGHT  FLOW. 

The  experimental  data  taken  during  testing  of  the  third  TF33  will  be 
evaluated  for  comparison  with  clean  air  baseline  data.  The  baseline  weight  flow  versus 
rotational  speed  plot  was  generated  from  the  first  three  clean  air  calibration  (run  nos. 
94,  95,  96).  A  fourth  order  polynomial  fit,  Figure  19,  was  performed  on  the  baseline 
data  with  a  99%  confidence  interval  of  a|w  /o~  /  6  I  =  4.38  lbs/sec.  The  deviation 
from  the  baseline  curve  was  then  evaluated  as  dust  ingestion  progresses.  It  is  expected 
that  as  ingestion  of  dust  progresses  the  weight  flow  of  air  through  the  engine  will 
decrease  with  respect  to  rotational  speed.  The  deviation  of  experimental  data  from 
the  baseline  curve  is  normalized  by  the  Z  statistic  corresponding  to  a  99%  confidence 
interval  for  an  assumed  standard  normal  distribution.  Figure  20  illustrates  the  normalized 
distribution  of  weight  flow  in  order  of  increasing  cumulative  dust  ingested.  A  significant 
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number  of  data  points  fall  outside  the  undeteriorated  99%  confidence  interval  indicating 
a  change  in  engine  performance  outside  the  computed  fluctuations  of  the  undeteriorated 
data. 

The  predicted  decrease  in  corrected  weight  flow  at  N2//o^  of  69.9%  speed 
is  at  a  normalized  deviation  of  1.71.  This  value  represents  the  average  expected  for 
a  sample  of  data  at  69.9%  rotational  speed.  By  assuming  that  the  standard  deviation 
of  the  data  for  the  deteriorated  engine  is  equal  to  the  standard  deviation  of  the 
undeteriorated  engine,  we  may  compute  the  99%  confidence  interval  of  the  deteriorated 
performance.  Note  that  the  predicted  value  corresponds  to  the  last  run  since  the 
measured  erosions  were  taken  from  the  blades  of  a  torn  down  engine.  Figure  21  shows 
the  predicted  deteriorated  99%  confidence  interval  and  associated  sampled  data  points 
for  the  last  few  runs.  The  majority  of  the  sampled  data  does  fall  within  the  99% 
confidence  interval  predicted  and  indicates  good  agreement  with  predicted  results. 
Some  of  the  dispersion  of  the  data  may  be  due  to  the  fact  that  the  majority  of  the 
data  were  sampled  at  the  design  speed  and  the  calculations  were  performed  for  the 
69.9%  speed  value.  It  is  not  possible  to  determine  the  functional  dependence  of 
deviation  from  the  baseline  curve  on  rotational  speed  due  to  the  limited  range  of  power 
settings  sampled. 
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SECTION  7 
CONCLUSION 


A  measurements  program  currently  underway  at  Arvin/Calspan  has  been 
used  in  the  evaluation  of  observed  engine  behavior  during  dust  ingestion.  The  major 
effect  of  solid  particle  ingestion  was  the  erosion  and  subsequent  performance 
deterioration  of  the  HPC.  The  most  noticeable  feature  of  blade  removal  was  increased 
tip  clearance  and  blade  leading  edge  erosion.  The  increase  in  proportionate  blade  loss 
through  the  compressor  will  result  in  substantial  performance  deterioration  in  both  the 
LPC  and  HPC  stages. 

The  observations  obtained  from  the  analysis  and  data  described  in  this 
report  indicate  that  a  deteriorated  engine  will  be  increasingly  likely  to  surge  when  low 
speed  acceleration  is  initiated.  Indications  of  deterioration  were  found  to  be  reduced 
burner  pressure  and  air  weight  flow  for  a  given  HPC  speed.  In  an  eroded  engine  before 
the  final  stage  of  deterioration,  the  intercompressor  bleed  and  anti-ice  air  served  to 
decrease  the  susceptibility  to  surge  at  low  power  settings.  Specifically,  use  of  anti-ice 
air  can  serve  to  decrease  the  likelihood  of  off-idle  surges;  use  of  intercompressor 
override  can  serve  to  '->id  surge  at  bleed  door  closing.  Off-idle  surges  emanate  from 
the  HPC  and  intercompressor  bleed  door  closing  surge  emanates  from  the  LPC.  The 
influence  of  erosion  is  to  make  the  front  stages  of  the  HPC  and  LPC  more  prone  to 
surge.  It  was  also  observed  that  near  cruise  power  the  fuel  controller  can  send  a 
deteriorated  engine,  while  operating  at  steady  state  and  having  a  narrow  surge  margin, 
into  a  autoacceleration  and  surge. 

To  quantify  the  performance  deterioration  a  predictive  capability  was 
developed  based  on  a  semiempirical  method.  The  individual  stages  were  modeled  to 
evaluate  the  effects  of  increased  tip  clearance  and  leading  edge  effects.  The  change 
in  total  pressure  rise  and  efficiency  for  the  undivided  stages  were  calculated  and  used 
to  develop  new  performance  maps.  The  stages  were  then  stacked  to  arrive  at  a 
deteriorated  compressor  map.  The  qualitative  effect  of  shifting  the  constant  speed 
lines  toward  reduced  mass  flow  and  lower  pressure  ratio  is  in  agreement  with 
experimental  observations.  The  quantitative  results  of  reduction  in  corrected  weight 
flow  gave  good  agreement  with  experimental  data.  The  reduction  in  surge  margin  at 
the  69.9%  N2  surge  line  was  predicted  to  be  51%.  Although  quantitative  measurements 


of  surge  margin  are  not  available,  the  observed  engine  behavior  does  agree  with  the 
expected  results  of  such  a  loss. 
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APPENDIX  A 

DEVELOPMENT  OF  ERODED  LIFT  COEFFICIENT 


For  contributions  due  to  camber,  the  classical  thin  airfoil  assumptions 
yield  equation  A—  1; 


where 


( cos  0-l)d0 


-  =  1/2  ( 1 -cos 0 ) 

c 


(A-l) 


—  =  slope  of  mean  camber  line 
dx 

By  examining  the  airfoil  section  as  it  erodes  at  the  leading  edge,  an  eroded  camber 
line  may  be  formulated  in  terms  of  the  clean  camber  line  and  a  change  in  chord  (see 
Figure  22). 
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(A-4) 


Since  the  inlet  fluid  angle  is  not  changed,  the  effect  of  the  increased  angle  of  the 
chord  line  is  felt  in  the  camber  angle  used  in  the  derivation  of  thin  airfoil  theory  (see 
Figure  23).  The  camber  line  used  is  defined  by  a  chord  line  in  the  horizontal  plane. 
As  the  chord  changes  the  camber  decreases  and  the  integration  is  taken  from  to 

rr  ,  where  9^  =  arc  cos  (1-2  A  c/c).  Equation  A-5  formulates  the  contribution  to 
lift  coefficient  from  the  eroded  camber  line. 
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To  compute  the  change  in  blade  lift  coefficient,  equation  A-5  is  subtracted 
from  equation  A-l  which  will  eliminate  biasing  felt  by  the  classical  thin  airfoil  theory 
assumptions  resulting  in  equation  A-6. 
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Only  small  changes  in  chord  are  of  interest  here  and  thus  the  above 
equations  can  be  simplified. 
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Equation  A-8  can  be  integrated  noting  that  E  (AC)  is  a  constant  for  any  given  eroded 
state  to  obtain  equation  A-9. 
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APPENDIX  B 

AN  APPROXIMATE  POTENTIAL  FLOW  SOLUTION  FOR  A  THIN 
HIGHLY  CAMBERED  AIRFOIL 


The  method  described  below  is  an  extension  to  classical  thin  airfoil  theory 
without  the  assumptions  of  small  angle  of  attack  and  small  camber  angle.  For  purposes 
of  evaluation  of  the  perturbation  velocity  (u,  w)  the  sheet  vortex  is  taken  to  lie  on  the 
x-axis,  in  the  range  0  <  x<c  ,  rather  than  on  the  line  z  =  z  (x).  This  is  the  only 
assumption  of  the  solution  and  makes  use  of  the  fact  that  ]  z(x)]«  c. 


The  boundary  condition  at  the  surface  of  the  camber  line  is  a  zero  normal  velocity 


v  +  w'  (S)  =  0 
«>,n 


(B- 1) 


From  Figure  25  equation  B-2  is  formulated 
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From  the  derivation  of  W(X)  (33)  we  have  equation  B-4. 
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(B-5) 


A  solution  for  y(A)  is  obtained  subject  to  the  Kutta  condition  y ( it')  =  0. 
This  solution  will  make  the  camber  line  a  streamline  of  the  flow.  The  solution  is 
presented  below  (33). 
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The  following  is  obtained 
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Note  the  similarity  to  the  Fourier  cosine  series  given  below. 
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where  f(0)  is  analogous  to 


where 


It  may  also  be  shown  that  the  lift  coefficient  of  the  airfoil  is  given  by 
equation  B-6  (33). 
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Substitute  into  equation  B-6. 


expanding  F  (x)and  noting 
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With  the  above  expression  and  a  knowledge  of  the  airfoil  camber  line  the  integral  can 
be  computed  numerically  and  the  approximate  potential  flow  solution  for  a  thin  highly 
cambered  airfoil  obtained. 


The  P/W  J57  15th  stage  blade  profile  was  fitted  to  the  following  polynomials. 

0  <  -  <  .13559 
“  c  ~~ 

-  =  -3.69107  fx/c)2  +  1.62518  (x/c) 
c 

.13559  <  x/c  <  .39831 

-  =  -1.37997  ( x/c )2  +  1.12402  fx/c)  +  .025466 
c 

.  398.31  <  x/c  1.0 

=  -.65337  fx/c)2  +  .49107  fx/c)  +  .16230 


z 

c 


The  solution  obtained  is  given  oelow 
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Figure  24.  Placement  of  the  vortex  sheet  for  thin  airfoil  analysis. 


Figure  25.  Determination  of  the  component  of  freestream  velocity  normal  to  the  camber  line. 
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APPENDIX  C 

DERIVATION  OF  COMPRESSOR  EXIT  PROPERTIES 


In  order  to  obtain  a  third  equation  and  thus  allow  the  calculation  of 
density,  velocity,  and  pressure  at  the  exit,  the  second  law  of  thermodynamics  and  the 
definition  of  polytropic  compression  efficiency  were  used. 
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Polytropic  compression  efficiency  enables  the  formulation  of  total  temperature  rise 
from  a  known  total  pressure  rise  and  overall  compressor  efficiency. 
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The  second  law  w,  ill  relate  the  local  static  to  local  stagnation  properties. 
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The  thermodynamic  equation  of  state  for  a  perfect  gas  will  relate  static  temperature 
to  static  pressure  and  density  and  allow  the  formulation  of  pressure  in  terms  of  density 
rise,  polytropic  efficiency,  and  total  pressure  rise. 


T  =  P/pR 


(C-6) 


Substituting  equations  C-k  to  C-6  into  C-2  and  rearranging  yields  equation  C-7. 
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V3 

Substituting  equation  C-7  and  continuity  into  the  definition  of  total  pressure  we  arrive 
at  an  expression  for  the  mean  axial  velocity  at  the  compressor  exit,  equation  C-8. 
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Where 


(T-Ofr-v) 


From  the  solution  of  velocity  the  density  may  be  computed  from  the  continuity  equation. 
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APPENDIX  D 
NOMENCLATURE 


B 

b 

CD 

C£ 

cL 

Cd 

Cx 

c 

DP 

EPR 

EGT 

E(  A c) 

f  (  ap) 

f 

g 

H 

HPC 

h 

K 

L 

LPT 


LPC 

MRT 

m 


erosion  velocity  exponent 

ratio  of  bleed  weight  flow  at  compressor  inlet 

particle  size  exponent 

cascade  drag  coefficient 

blade  lift  coefficient 

cascade  lift  coefficient 

particle  drag  coefficient 

axial  velocity 

chord 

particle  diameter 

engine  pressure  ratio 

exhaust  gas  temperature 

leading  edge  coefficient 

dust  angle  of  approach  function 

fuel-air  ratio  at  combustor  inlet 

erosion  function 

total  or  stagnation  enthalpy,  Btu/lb 
high  pressure  compressor 
blade  height 

cascade  modifying  coefficient 
cumulative  mass  of  dust  ingested 
low  pressure  turbine 
low  pressure  compressor 
military  rated  thrust 
mass  flow,  (slugs/sec) 
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NOMENCLATURE  (Continued) 


Ni 

N2 

N 

P 

Pt 

PR 

PT7 

PLA 

PRS 

T 

rt 

Re 

eP 

rP 

R 

SM 

s 

T 

TRT 

Tt 

u 

U 

v 

W 

w 

wm 

x 


percent  of  LPC  rotational  speed 

percent  of  HPC  rotational  speed 

number  of  stages 

static  pressure 

total  pressure 

total  pressure  ratio 

total  pressure  at  turbine  exit 

power  level  angle 

surge  pressure  ratio 

hub-tip  mean  radius 

tip  radius 

particle  Reynolds  number 
mean  radius  of  particle 
gas  constant 
surge  margin 
cascade  blade  spacing 
static  temperature 
takeoff  rated  thrust 
total  temperature 
velocity 

rotational  speed 
velocity 

vertical  velocity  at  camber  line 
weight  flow 

mean  relative  fluid  velocity 
distance  along  chord  line 


68 


NOMENCLATURE  (Continued) 


height  of  camber  line 
airfoil  angle  of  attack 
relative  fluid  angle 
ratio  of  specific  heats 

ratio  of  total  pressure  to  NASA  standard  sea-level  pressure  at  2116  psi 

volume  removed  per  particle  impact 

efficiency 

ratio  of  total  temperature  to  NASA  standard  sea  level  temperature  of 
518.7°R 

tip  clearance/blade  height 

viscosity 

stagger  angle 

total  pressure  loss  coefficient  wrt.  tip  speed 
density 

reduced  density 

power  of  fan  and/or  accessories,  Btu/s 

tip  clearance 

particle  relaxation  time 

flow  coefficient 

total  pressure  coefficient 

fraction  of  blade  height  influenced 

Subscripts 
before  impact 
after  impact 
compressor  inlet 
compressor  outlet 
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NOMENCLATURE  (Continued) 


i 

ii 
is 
L 
m 
N 
o 

P 

pc 

T 

t 


cascade  inlet 
cascade  outlet 
isentropic 
loss 

cascade  mean  fluid  value 

norma! 

uneroded 

particle 

polytropic 

tangential 

tip 

free  stream 
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